Surface Code Threshold in the Presence of Correlated Errors 
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We study the fidelity of the surface code in the presence of correlated errors induced by the 
coupling of physical qubits to a bosonic environment. By mapping the time evolution of the system 
after one quantum error correction cycle onto a statistical spin model, we show that the existence of 
an error threshold is related to the appearance of an order-disorder phase transition in the statistical 
model in the thermodynamic limit. This allows us to relate the error threshold to bath parameters 
and to the spatial range of the correlated errors. 



The surface code is considered one of the best quan- 
tum error correction codes (QEC) to implement on phys- 
ical devices This steams from two major points: 
first, all syndromes and operations can be performed with 
spatially local operators; second, all threshold estimates 
show that, for sufficiently large lattices, the error thresh- 
old is the highest known for two-dimensional architec- 
tures with only nearest- neighbor interactions [|], |[. 

The error threshold is usually defined for stochastic 
error models. By assuming that errors are independent 
events and assigning a probability p to each of these 
events, it has been shown numerically that the quantum 
information encoded in the surface can be faithfully pro- 
tected when p is below a critical value. Although this 
result is firmly established numerically, there are two big 
open questions that still need to be addressed. First, 
stochastic error models are approximations to reality that 
sometimes cannot be justified. In fact, most studies so 
far lacked a microscopic description of the interaction be- 
tween physical qubits and the environment. Second, the 
same locality of operations and syndromes that makes 
the surface code powerful also makes it more susceptible 
to correlated errors. Thus, a discussion of the tradeoff 
between locality of operations and correlated errors is 
long overdue. In this Letter, we address both issues by 
employing a more realistic error model. We consider a 
Caldeira-Leggett type of environment where freely prop- 
agating bosonic modes couple linearly and locally to the 
physical qubits. Such a model has a very strong physical 
motivation since in most experimental implementations 
photons and phonons couple to the two-level systems 
making up the physical qubits. This model also plays 
a fundamental role in our understanding of decoherence 
§ and its interplay with QEC f^Q. 

Consider a logical qubit in a quantum memory. In a 
QEC cycle, the logical qubit is prepared and left to freely 
evolve during a certain time interval. Then, syndromes 
are extracted and, if necessary, suitable error correction 
operations are implemented to bring the logical qubit 
back to its original state. In our analysis, we evaluate the 
fidelity of a logical qubit after such a QEC cycle under 
the assumption of nonerror syndromes. This assumption 



is not essential to our results, but makes the calculation 
more concise. In addition, we assume that the bath is 
initially at zero temperature and that it is reset to this 
temperature at the end of the QEC cycle. 

Our results show a sharp transition between two dis- 
tinct noise regimes. On one hand, below a fictitious crit- 
ical "temperature" (which is related to microscopic pa- 
rameters of the bath), the dissipation due to the bath 
cannot be suppressed by the encoding. On the other 
hand, above this critical "temperature", we show that 
if the thermodynamic limit is taken, the effects of the 
bosonic environment become irrelevant and the logical 
qubit is fully protected. 

Even though this investigation focuses on quantum in- 
formation protection, the physical problem we consider 
has a much broader appeal. In essence, it amounts to 
a lattice gauge system interacting with a scalar bosonic 
field in two dimensions |l^] . In this language, the nonex- 
istence of a quantum error threshold can be understood 
as the lifting of the ground state topological degener- 
acy. Our discussion of the error threshold can therefore 
be recast as a quantum phase transition. This fits into 
our earlier discussion of the threshold theorem as resem- 
bling a quantum phase transition [pi. This analogy is 
nontrivial since the error threshold is in essence a driven 
dynamical problem, far from the equilibrium conditions 
required to describe phase transitions in statistical me- 
chanics. This work turns what was an analogy into a well 
defined map. Therefore, we believe that the results pre- 
sented here transcend our original motivation and com- 
plement the recent discussion in Ref. 

The model - We consider physical qubits {oi}-j=;L,...,iv 
(i.e., spin 1/2 systems) located on the links of a square 
lattice. The QEC code is defined by an encoding pre- 
scription and a set of stabilizer operators The sta- 
bilizer operators of the surface code are easily labeled 
when we define stars and plaquette operators (see Fig. 
|l|). Star operators are the product of the four a x op- 
erators of qubits adjacent to a vertex of the square lat- 
tice, Ao — Yiie<> a i- Similarly, plaquette operators are 
the product of four a z of quits located at the edges of 
a tile of the square lattice, £?□ = ILea^f ■ Each pla- 
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FIG. 1. (a) Star (A) and plaquette (_B) operators, (b) Ex- 
amples of logical X and Z operators. The crosses and shades 
indicate the qubits involved in the respective operator. 



quette and star is a stabilizer that has to be measured 
at the end of a QEC period. The logical operators are 
defined as X = Jlier^f' where T is a path along the 
center of the plaquettes crossing horizontally the lattice, 
and Z = Jligf a ii wnere f is a path along the edge of 
the plaquettes crossing the lattice vertically. Finally, the 
codewords of the surface code are also easily written, 

{\\)=G\F Z )M)=X\\)}, (1) 

Y[a, (1 + A^), Na is the total number 



where G = , ■ . i 
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of possible stars, and \F Z ) is the ferromagnetic state along 
the positive z direction, namely, \F Z ) = Yl i=1 | t)i,«> 
The Hamiltonian we consider is written as 



H = H + V, 



(2) 



where H is a free bosonic Hamiltonian, H a — 
D k -£ Wka^ak, and V = f 2 r /(r)of, where r denotes 
the spatial location of a qubit and / is a local bosonic 
operator, 
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Here, D is the bath spatial dimension, w k = u|k|, and ui 
is a microscopic coupling that makes / dimensionless (we 
adopt units such that h = 1). This error model leads to a 
remarkably simple evolution operator in the interaction 
picture: for a time interval A, we have JllJ 



f/(A) = e 



- p-V we rr (A) 



' £^ 8 *rs(A)« 
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(4) 



where N is the total number of qubits, Q, 
(0|F r (A)F s (A)|0), F r (A) = f n A dt f( T ,t), $ 



(A) = 

f*dtf(r,t), $ rs (A) = 

G rs (A)+ J Q a dti Jq 1 dt 2 [f(r,t 1 ),f(s,t 2 )}, and :: stands for 
normal ordering. For ohmic baths, the correlation func- 
tion for r/s takes the simple form [pi 
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f , < |r-s| < vA, 
, < vA < |r - s|. 

(5) 



Thus, we can introduce a fictitious inverse tempera- 
ture P = 5 (^ !a ) and rewrite the intermediate factor 
that contains the two-spin interaction in Eq. (Eh as 
e p ^r^s J ™ r s ; w here J rs represents an effective anti- 
ferromagnetic coupling. 

To simplify the notation, we assume that the system 
is prepared initially in the logical state |f) and the boson 
field initial state is the vacuum, 



IVto) = (G\F Z )) ® |0). 



(6) 



We then let the system evolve under the Hamiltonian 
H until a time A, when an error correction protocol is 
performed flawlessly. 

The assumptions - Since we are only allowing for bit- 
flip errors, the syndrome outcome for the star stabilizers 
is trivial. For the plaquettes, in principle, all possible 
syndromes should be considered. However, it is useful 
to look at the most benign evolution and assume that 
all plaquette syndromes return a nonerror. This non- 
error syndrome provides an upper bound to the available 
computational time and also substantially simplify the 
calculation since it removes from consideration which re- 
covery operation should be performed to steer the system 
back to the computational basis. 

In QEC theory, it is standard to focus only on the 
qubit system evolution and disregard any change to the 
environment's state, even though the latter is a also a 
quantum system capable of sustaining correlations. If no 
extra step is taken to dissipate those correlations, the en- 
vironment will keep a memory of events that happened 
between the QEC periods. Keeping track of such bath- 
induced, long-time correlations between QEC cycles in 
a fidelity calculation is a difficult task even for simple, 
nontopological logical qubit systems 0, O. Thus, to 
proceed with the calculation, we consider an extra step 
to the QEC protocol. In addition to projecting the quan- 
tum computer wave function back to the logical Hilbert 
space, we assume that at the end of the QEC step the 
environment is reset to its ground state. Hence, from 
this point on, we are excluding from the calculation any 
spatial correlation between QEC periods, as well as mem- 
ory and spatial correlations between the time evolution 
of bras and kets. When we adopt this extra simplify- 
ing assumption, we can conveniently rewrite the nonerror 
syndrome projector as 



P' = \iPo)(i>o\+X\^ )(i;o\X. 



(7) 



Using Eq. (Q), it becomes now straightforward to write 
an expression for the logical qubit fidelity just after the 
syndrome extraction, 



T = 



\A\ 



where A = (ip \U(A)\ip } and B = (ip \XU(A)\ip Q ) 



(8) 
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Fidelity calculation - Thus, our task now is reduced to 
evaluate A and B. Using Eq. (^), it is straightforward to 
show that 



and 



where 



A = X (F z \e-^G 2 \F z ), 



B = X (F z \Xe-^\F z ), 



(9) 



(10) 



(11) 



and x — e~ 2 r A ' err ( A ). Notice that although we chose 
to start with a microscopic model of the environment to 
make a connection to physical implementations, we could 
as well have started by imposing an effective two-body 
interaction between qubits such as that defined by %. 

Clearly, when (3 — > 0, we have B — > and J- — > 1. A 
perturbative expansion for small /3 is the standard route 
to discuss the error threshold. This "high-temperature" 
expansion will be discussed elsewhere. Here we follow a 
different route. To understand the opposite limit, /3 — ¥ 
oo, we need to rewrite \F Z ) in the x basis, 
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This dual representation corresponds to a "low- 
temperature" expansion and it is suitable for describing 
the regime where error correlations are strong. Inserting 
Eq. @ into (g) and @, we obtain 

X 
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- f3E °(S\G 2 \S), 



(13) 



(S\XG 2 \S), 



(14) 



respectively, where 15) is an element of the x basis and 
E s = (S\H\S) is its "energy". Notice that E 8 may have 
an imaginary part. 

If we had unrestricted sums in Eqs. (|lf) and @, we 
would be essentially discussing a two-dimensional Ising 
model. However G 2 projects 15*) onto the subspace of 
positive stars, and among other things remove time re- 
versal state of \S±) from the sum. This restriction makes 
the computation of A and B nonstandard. Nevertheless, 
the action of X is to introduce a sign between two dis- 
tinct classes of states. To better understand this, we need 
a more convenient way to write the states \S) in the re- 
stricted subspace of positive stars. It is not hard to prove 
that these states fall into two groups, and 
where 



\S 4 



U.B a .\F x 



and 



\S-) 



(15) 



(See Fig. ||.) Here, JTj B\j. is a product of plaque- 
ttes that do not touch a logical error Z 1 and \F X ) is 
the ferromagnetic state in the x basis. Splitting the 
terms of the sums in Eqs. (13) and (|l4| ) between these 



two groups of states, we rewrite A = ^w(T+ + 71) and 
B= 2^(T+ -T-), where 



7i = E<^ 

s± 



' PH \S±) 



(16) 



Thus, we have to evaluate the "free energies" of the two 
groups of states (see Fig. ||) . 

The double sum in % [see Eq. (pi])] runs over lat- 
tice points r and s inside and outside the path of the 
logical operator Z 1 . Thus, let us break these points 
into two sets, namely, {r} = {t 7 } {u 7 } and {s} = 
{v 7 } {w 7 }, where t 7 and v 7 belong to Z 7 while u 7 
and w 7 do not. Furthermore, let us factor the sum over 
all S-\- into a sum over paths 7 and a sum over con- 
figurations Sj_ compatible with a logical operator along 
this path, namely, T^s± = S 7 ^2s]_ • After some sim- 



ple algebra, we obtain 7± = X) 7 
^ = {Sl\Et v Jtvof<|^> and 
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where 



(17) 



with hi, = J2t ^ tw - ^ e can see ^ na ^ t ne effect of the 
logical operator Z 7 is to introduce a boundary term rep- 
resented by the effective local magnetic field h^,. 

We are now in the position to state our definition of 
the quantum error threshold. We define the critical pa- 
rameter flc as the value of j3 that separates the regime 
where J- = 1 from the regime where J- < 1 in the ther- 
modynamic limit (N — > 00). 

Phase transition - The evaluation of z± in Eq. (^) is a 
formidable task and a general answer may only be achiev- 
able through numerical simulations. Hence, we now re- 
strict our considerations to more manageable effective 
qubit interactions. Let us first consider the case of 



J, r, s nearest neighbors, 
otherwise, 



(18) 



where J is real. Such case is of physical relevance. Any 
measurement or operation on the plaqucttcs and stars 
can introduce short-range correlated errors, regardless of 
the presence of an environmental bath. For an ohmic 
bath, it corresponds to set vA to the order of the lattice 
spacing. 

To show that an order-to-disorder transition indeed ex- 
ists, let us consider Eq. (17) with (|l8| ) in the absence of 
the field h%,. Using Eq. (15), we obtain 
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FIG. 2. Example of states that contribute to the S+ (a) and 
S- (b) sums [see Eq. ([[5])]. The light (dark) circles indicate 
qubits in + ( — ) a x eigenstates. The dark line indicates the 
path of the logical error , which is active in (b) and inactive 
in (a). 



where we have multiplied all spins in one of the sublat- 
tices by —1 to make the model ferromagnetic. The prod- 
ucts over plaquettes introduce a sum over loops where 
the qubits at the these loops have negative eigenvalues 
(they are located at the edges of the shaded regions in 
Fig. ||). Let us label these qubits by a and b and use c 
and d to label qubits with positive eigenvalues. Thus 



Jed P ~I3 Ea 



(20) 



At "low temperatures" , the first two factors in Eq. ( pQ) 
dominate. The third factor leads to excitations above 
the ferromagnetic ground state whose "energy" is Ep. 
Therefore, we can go back Eq. ([l7]) and write (for — 
0) 



-2/3 E, 



£loop,c^loop 



(21) 



loops 



Then, the "energy" cost of a loop with length £ is equal 
to 4 J£. The scaling of the number of loops with £ is well 
known: Ng ~ fjr , where fi is the lattice connectivity con- 
stant (e.g., fx ~ 2.638 for a square lattice) Hence, 
the contribution of loops with length I to J2sj_ zl § oes 
as // e~ 8l3je . Clearly, a tradeoff between "energy" mini- 
mization and "entropy" maximization takes place and a 
phase transition is expected at f3 = /3 C , where 



In fi 
~8J' 



(22) 



For f3 > (3 C the system is in its ordered phase. Thus, we 
expect the system to be sensible to the boundary mag- 
netic field direction, leading to T < 1 after the QEC 
period. In particular, there is an obvious result to be 
stated. At the zero "temperature" limit, f3 — > oo, only 
the "ground state" contributes to the sums (i.e. only the 
T- sum survive). Hence, = \B\ and T = |. Con- 
versely, for (3 < f3 c (therefore A < A c ), a disordered phase 



develops, leading to the insensibility to the boundary field 
and, therefore, T = 1. 

There are some roadblocks to study the general case 
represented by the correlator in Eq. (|j). The first one 
is its imaginary part. It is expected that any imaginary 
part would have a negative effect on the fidelity (a conclu- 
sion that can be reached by applying Schwarz inequality) . 
Even if we disregard this imaginary part, the magnetic 
model given by the real part of Eq. (0) may lead to frus- 
tration. However, since the real part of the correlator 
is a slow growing function, it is likely that the model is 
controlled by boundary effects (as many other features of 
topological systems are). Only through future numerical 
work it will be possible do address these issues. 

In order to gain some insight into the effects of long- 
range correlations, we arbitrate a stripped version of the 
ohmic model. Consider the case where qubits on one of 
the sublattices interacts only with qubits on the other 
sublattice thought the real part of Eq. (|^) . This interac- 
tion preserves the bipartite nature of the lattice and we 
can again map it onto a ferromagnetic model. We can 
repeat all step of the calculation done in the case with 
nearest neighbor interactions. Then, the "energy" cost 
of a loop with length I is equal to (n — 2)J£, where n is 
the number of sites connected to any lattice site through 
the interaction (i.e., n is related to the range of the in- 
teraction). The transition happens at 



In /i 



2(n-2)J 



(23) 



Note that j3 c explicitly depends on the interaction range. 
In a strict sense, if all qubits in the lattice are in the 
"causality cone" (and therefore participate in H) there 
is never a threshold. Although this could seem a dismal 
result, it also shows us that a a finite QEC period intro- 
duces an infrared cutoff that creates a finite transition 
"temperature" . The smallest of such transition "tem- 
peratures" corresponds to the case of nearest neighbors, 
Eq. (|l8|). Finally, we note that, under the same assump- 
tions that we used for the ohmic bath, superohmic envi- 
ronments are local and therefore will always yield a finite 
transition "temperature" . Conversely, subohmic baths 
will require an infrared cutoff to yield a finite threshold. 

In conclusion, in this Letter we studied the fidelity of 
a logical qubit encoded with the surface code after one 
complete QEC period. We derived a nontrivial mapping 
to an statistical mechanical problem and provided an an- 
alytical expression for the error threshold for some corre- 
lated models. This mapping provides a promising route 
for exploring fault tolerance of topological quantum error 
correction codes in the presence of realistic environments. 
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